arXiv: 1509.0501 Ivl [astro-ph.SR] 16 Sep 2015 


Report on the O IV and S IV lines observed by IRIS 

Peter R. Young 

pyoung9@gmu.edu 

George Mason University & NASA Goddard Space Flight Genter 


This report has not been submitted to a journal but is provided on astro-ph as a reference for 
other scientists who use IRIS data. 

The Oiv intercombination lines observed by IRIS between 1397 and 1407 A provide useful 
density diagnostics. This document presents data that address two issues related to these lines: 

1. the contribution of Siv to the Oiv A1404.8 line; and 

2. the range of sensitivity of the Oiv A1399.8/A1401.2 ratio. 

1. Overview 

Oiv has a ground conhguration of 2s^2p, and Siv a ground conhguration of 3s^3p and this 
similarity leads to common features in their spectra. In particular, each ion has a set of in¬ 
tercombination transitions of the form ^P-^P where the upper state lies in the excited nsnp^ 
conhguration. There are hve distinct transitions in these multiplets, and the O iv lines he be¬ 
tween 1397 and 1408 A, and the Siv transitions he between 1398 and 1424 A (Tabled]). The 
O IV multiple! is signihcantly stronger than the S iv multiple! in most conditions, and both sets 
of multiplets show density sensitivity in the log iVe=10-13 range. 

Table 1: Oiv and Siv transitions. 


Ion 

Wavelength 

Transition 

Gomment 

Oiv 

1397.198 

1/2-3/2 

weak line 


1399.766 

1/2-1/2 



1401.157 

3/2-5/2 



1404.779 

3/2-3/2 

blended with S iv 


1407.374 

3/2-1/2 

not observed by IRIS 

Siv 

1398.040 

1/2-3/2 

very weak line 


1404.808 

1/2-1/2 

blended with 0 iv 


1406.016 

3/2-5/2 



1416.887 

3/2-3/2 

not observed by IRIS 


1423.839 

1/2-3/2 

not observed by IRIS 


Since S iv is formed very close in temperature to O iv, then generally there is little interest 
in studying the S iv lines over the stronger O iv lines. However, one of the S iv lines blends with 
a useful Oiv transitio n and so it is nece s sary t o take account of this (Sect. [5|). 

The IRIS satellite ( De Pontieu et ah . 2014 1 observes the 1389.0-1407.0 A wavelength band, 
although in practice the upper limit is ~ 1406.6 A. This means that a number of the transitions 
are not observed by IRIS (see Tabled) for details). 

There have been many previous studies of the O iv and S iv lines due to their importance in 
both solar and stellar spectra, and Table [2] lists some key papers. 
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Table 2: Previous studies of the O IV intercombination lines. 


Paper 

Cook et al. 119951 


Harper et al. (19991 
Teriaca et al. (2001 1 
Keenan et al. (2002 1 


Keenan et al. 


Comment 

Oiv and Siv lines in HRTS solar spectra and HST/GHRS 
spectra of Capella. 

Oiv and Siv lines in HST/GHRS spectra of RR Tel. 
Densities in explosive events from SOHO/SUMER. 

Oiv and Siv lines in HST/STIS spectra of RR Tel, and 
SOHO/SUMER solar spectra. 

OIV lines in stellar and solar spectra. 


Atomic data 


Atomic data for Oiv and S i v used in this work are obta ined from version 8 of the CHIANTI 
atomic database ( Dere et ah . 1997 : Del Zanna et ah . 2015ll . Radiative decay rates for the lowest 
20 atomic lev els (including the levels g iving rise to the intercombination lines) of O iv were 
obtain ed from Correge fc Hibbert ( 20041 ). and all other decay rates were ta ken from Liang et al.l 


although an error was corrected in these data, as reported in iDel Zanna et al 


( 2 OI 2 I I. Effective collision strengths for all transitions wer e taken from Liang et al. ~ |2012l l 


The Siv data-set has been unchanged since CH I ANTI 5 dhandiet al.l, 


200 ' 


(|201,5l l. 

and consists 


of radiative decay rates from Hibbert et al. (2002), Taval ( 1999ll . I Johnson et al. (1986) and 


some unpublished d ata of P.R. Young. The data for the intercomb ination tran sitions are from 
Hibbert et al. (2002). Effective collision strengths for Siv are from Taval ( 2 OOC 1 II . 


3. O IV Density sensitivity 

Figured] shows the density sensitivity of the three Oiv density diagnostic ratios that are impor¬ 
tant for IRIS. Temperature sensitivity is relatively small close to the temperature of maximum 
ionization of Oiv (log T = 5.15), but in Figured) we show the ratios at log T = 4.40 at which 
there are relatively large differences. It is possible in non-equil i briuin conditions that O iv could 
be formed closer to chromospheric temperatures ( Olluri et ah . 2013l l. 


4. IRIS spectra and calibration 

In order to derive densities from the Oiv density diagnostics, a number of spectra were selected 
and these are identified in Table [3l The spectra were chosen based on criteria such as: strength 
of the lines; lack of blending; symmetry of lines; and possibility of high density. The spectra 
were created using the IRIS_SUM_SPEC routine in Solarsoft. This routine sums blocks of pixels 
to create averaged spectra with error bars. For example, the IDL call to generate spectrum no. 8 
is: 


IDL> file=iris_find_file(02-oct-2013 21:34’) 

IDL> str=iris_sum_spec(file,xpix=9,ypix= [113,115],rnum=l) 

The X-pixel and Y-pixel values for each spectra are given in Table [S] as well as the raster 
number (for data-sets in which multiple rasters were obtained). For each spectrum a background 
spectrum was subtracted from the feature spectrum by identifying a suitable background location 
near the feature along the slit direction. The Y-pixel values for the background spectra are 
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Figure 1: Oiv ratio plots generated with CHIANTI 8. The black line has been calculated at 
log T = 5.15 and the red line at log T = 4.4. 

given in Table E] (the column “BG Ypix”). Note that the same X-pixel values as for the feature 
spectrum were used. 

The background spectrum was subtracted from the feature spectrum, with the error arrays 
added in quadrature. Emission lines in the subtracted spectra were then fit with Gaussians using 
the IDL routine SPEC_GAUSS_WIDGET, which is also available in Solarsoft. Line intensities 
in data number (DN) units derived from the Gaussian fits are given in Table 01 

In this work we simply use the line intensities given in DN. The current version of the 
IRIS instrument response (implemented through the IDL routine iris_get_response) gives a 
completely flat response through the wavelength range of the O iv lines, and so the response 
function has no impact on line ratios. Note that densities are computed by using theoretical 
ratios given in photon units rather than energy units. 


Table 3; Parameters used for computing spectra. 


Index 

File start time 

Rast. 

no. 

Xpix 

Ypix 

BG 

Ypix 

Comment 

1 

ll-Oct-2013 23:54 

0 

238 

471:475 

460:464 

SIV strong and fairly symmetric 

2 

ll-Oct-2013 23:54 

0 

76:79 

714:718 

733:737 

Supersonic loop footpoint; fit 
only supersonic component 

3 

ll-Oct-2013 23:54 

0 

64 

876:881 

914:919 

Bright, elongated structure 

4 

ll-Oct-2013 23:54 

0 

92 

710:714 

658:662 

Loop footpoint, not supersonic 

5 

lO-Sep-2014 11:28 

0 

2321 

338:341 

344:346 

Flare kernel site, narrow line 

6 

lO-Sep-2014 11:28 

0 

1988:1992 

491:492 

481:485 

Long loop structure; lines have 
weak, extended wings 

7 

lO-Sep-2014 11:28 

0 

2381 

301:304 

309:313 

Flare kernel, narrow lines, ex¬ 
tended red wing 

8 

24-May-2014 06:10 

0 

25:32 

532:545 

— 

South coronal hole, just above 
limb 
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The line intensities derived from the Gaussian fits are given in Table 01 and properties 
derived from the fits are given in Table O For the Oiv A1399.8/A1401.2 ratio, we give the 
measured ratio and the corresponding density, derived from CHIANTI 8 assuming a temperature 
of log T = 5.15. The measured wavelength sep aration of the two l ines is given in Table [5] which 
can be compared with the separation given bv lYonng et al.l ( 20111 ) of 1.391 A. We also give the 
ratio of the width of A1399.8 to that of A1401.2, which gives an indication of the accuracy of the 
fits since the two lines should have the same widths. Large differences may indicate unaccounted 
for blending. The final column in Table [5] is the percentage contribution of S iv to the line at 
1404.8 A and this is discussed in the following section. 


5. The 1404.8 line and S IV 

The Oiv A1404.8 line is potentially a good density diagnostic relative to A1401.1 (Fig. [T]), but 
there is a known blend with Siv A1404.8. The Siv contribution can be estimated through a 
branching ratio with A1423.8, but this line is not observed by IRIS. An estimate can be made 
by considering A1406.0, an d this is discussed be l ow. 

Based on Skylah data, Feldman fc Doschek ( 1979l i found that Siv contributed 10% to the 


observ ed feature in quiet Sun conditions, but this rises to 33% in flare conditions. ICook et al 


(Il99.5li state that S iv contribut es less than 1 0 % by comparing the strength of the line against 
the nearby A1406.1 transition. Teriaca et al. ( 200ll ) stated that Siv only contributed 3.6% to 
the observed feature in SUMER spectra but no details were given. 

As IRIS does not observe the S iv A142 3.8 line, then we only ha ve the A1406.0 line to estimate 


the Siv contribution. Using CHIANTI 8 (jPel Zanna et al.l . 1201.51 ). we find the A1404.8/A1406.0 


ratio increases slightly with density, with values of 0.199, 0.202, and 0.224 at log Ne=9, 10, and 
11, respectively. In Table[5]we multiply the A1406 intensity, where available, by 0.210 and then 
divide this quantity by the A1404.8 intensity to e stima te the percentage contribution of Siv to 
the observed feature. We note that I Cook et al.l (119951 ) used the same method and they found 
the same correction factor (0.21) despite using much older atomic data. 

The results demonstrate that Siv can make a significant contribution to the blended line 
and so we generally recommend not to use the 1404.8 A line as part of a density diagnostic, 
unless the S iv A1406.0 is also measured. We note, for example, that the IRIS flare line list does 
not include S IV A1406.0. The A1404.8 will be useful in low density data-sets, such as coronal 
holes and quiet Sun as A1399.8 is weak in these regions, and so the A1401.2/A1404.8 ratio should 
give more precise results. 


6. Blending lines for O IV 1399.8 

In this and the following section we consider blending lines for OIV A1399.8 and A1401.2 that 
may affect density measurements from the A1401.2/A1399.8 ratio. 

A1399.8 is partly blended with a Fell line at 1399.960 A and Fig. [2] shows an example where 
this line is quite strong, adding a narrow peak to the long wavelength side of the Oiv line. 
Visually the line is about the same strength as Fell A1401.777 which is close to Oiv A1401.2 
(see Figure [2]), and so this unblended line may be used to estimate whether the Fell contribution 
is important to Oiv. 

A more complex blending scenario can occur in flare ribbons, and Figure [3] shows an example 
spectrum for the 2014 September 10 X-flare in the vicinity of the Oiv A1399 line. The bright 
blob in the image is the O iv line, and one can see a number of very narrow features in the 
spectrum that have a larger extent in the Y-direction than the Oiv line. The strongest of the 
displayed lines is at 1399.69 A and, based on the intensity distribution along the slit, appears to 
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Table 4: Line intensities (DN). 




Oiv 


Siv 

Spectrum 

A1399 

A1401 

A1404'’ 

A1406 

1 

2576 ± 21 

7499 ± 29 

2239 ± 19 

2738 ± 24 

2 

600 ±8 

2820 ± 11 

1102 ±9 

179 ±5 

3 

511 ± 13 

1532 ± 15 

622 ± 14 

856 ± 14 

4 

256 ± 12 

1093 ± 14 

260 ± 13 

62 ±8 

5 

10375 ± 42 

25491 ± 60 

9750 ± 43 

— 

6 

184 ±6 

627 ±8 

148 ±6 

— 

7 

847 ± 15 

2119 ± 19 

594 ± 15 

— 

8 

142 ±3 

782 ±3 

360 ±3 

119 ±2 


Table 5: 

Derived properties. 


Oiv A1399.8/A1401.2 


Spectrum AA Intensity ratio Log Density Width ratio S iv % 


1 

1.379 

± 

0.001 

0.343 

± 

0.003 

1 1 99+0-02 
ll.ZZ_Q Q3 

1.030 

± 

0.009 

25.7 

2 

1.385 

± 

0.001 

0.213 

± 

0.003 

10 16+004 

ru.iD_g Q 3 

0.993 

± 

0.014 

3.4 

3 

1.386 

± 

0.002 

0.334 

± 

0.009 


0.994 

± 

0.027 

28.9 

4 

1.388 

± 

0.003 

0.234 

± 

0.012 

10.39l°;J° 

0.931 

± 

0.046 

5.0 

5 

1.387 

± 

0.000 

0.407 

± 

0.002 

12.02l°;06 

1.022 

± 

0.004 

— 

6 

1.389 

± 

0.002 

0.293 

± 

0.011 

10.85l°;0? 

0.897 

± 

0.032 

— 

7 

1.386 

± 

0.001 

0.400 

± 

0.008 

11-86+^3 

1.030 

± 

0.020 

— 

8 

1.385 

± 

0.001 

0.182 

± 

0.003 

9.55l°;J° 

0.944 

± 

0.017 

7.0 


be a H 2 line but the identification is not known. Care must be taken when deriving the intensity 
of Oiv A1399.8 from flare ribbons to estimate the strength of both this line and Fell A1399.96. 


7. Blending lines for OIV 1401.1 

A commonly-seen line close to A1401.1 is Si A1401.514, although it is well-separated from the 
Oiv line and generally not a problem (Fig. [2]). Note that the velocity of the Si line relative to 
the O IV line is 76.4 km s ~h Near sunspots strong supersonic downflows of ~ 90 km s“^ are 
often seen in the O iv line ( Straus et all . 20151 ). and so the S i line can potentially be a problem 
when studying such flows, however the author’s experience is that S i is generally negligible in 
these structures. 


8. The high density limit 

Figure [4] shows the densities derived from the Oiv A1399.8/A1401.2 ratio for each of the spectra 
identified in Table [3l The data span most of the density range of the ratio, from the coronal hole 
observation (spectrum 8) to the flare kernels (spectra 5 and 7). The fact that the measured ratios 
lie within the expected range of variation gives confidence in the atomic data for Oiv. The 2014 
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Figure 2: Spectra from an exposure obtained by a raster beginning at 22:08 UT on 2013 October 
12. The left two panels show detector images around the Oiv A1399.8 (upper) and A1401.2 
(lower) lines, and the middle and right panel show spectra obtained at Y-pixel 114 (indicated 
by the short blue lines on the images. 


September 10 X-flare does show exa mples of densities very close to 10^^ cm“^. We note that 
in features referred to as “bombs” bv iPeter et al.l (120141 ). and also in flare kernels the O iv lines 
often become very weak or disappear. This was noted from Skylab data bv iFeldman fc Doschek 
( 19781 ) and was suggested to be due to densities that are so high that the Oiv emitting levels 
are de-populated by electron collisions, reducing the strength of the lines. This would happen 
at ~ 10^^ cm“^. The fact that densities of ~ 10^^ can be measured with the Oiv ratio in later 
stages of flare kernel evolution suggests that higher densities may be feasible in the early stages 
of flares. 

As a sanity check on the extreme ratio values that were found in this analysis. Appendix [Al 
shows over-plots of the two O iv on top of each other. 


9. Summary 

The O IV density diagnostics observed by IRIS have been investigated and the following results 
found. 

1. The line observed at 1404.8 A is dominated by Oiv but can contain a contribution from 
S IV up to 29% and so it is recommended that the S iv A1406.0 line is observed in order to 
estimate the contribution. 

2. The OIV A1399.8 line is blended with an unknown line at 1399.69 A, which is likely due 
to H 2 . This line can be stronger than the OiV line in flare kernels. 

3. A low density of log Ng = 9.5 is measured just above the limb in a coronal hole. 

4. In two flare kernel sites of the 2014 September 10 X-flare the density reached log Ng = 12. 


References 

Cook, J. W., Keenan, F. P., Dufton, P. L., Kingston, A. E., Pradhan, A. K., Zhang, H. L., 
Doyle, J. G., & Hayes, M. A. 1995, ApJ, 444, 936 

Correge, G., &: Hibbert, A. 2004, Atomic Data and Nuclear Data Tables, 86, 19 

De Pontieu, B., et al. 2014, Sol. Phys., 289, 2733 


6 
























2500 

2000 

z 

Q 

-t 1500 

c 

B 

1000 


1399.5 1400.0 ^ 1400.5 

Wavelength / A 

Figure 3: Spectrum for the 2014 September 10 X-flare. An exposure from X-pixel 2314 is shown 
in the left panel, and has been saturated to better show the weak lines. The Oiv A1399 line is 
the broad, bright structure. The blue line indicates the ID spectrum shown in the right panel, 
which corresponds to Y-pixel 353. 

Del Zanna, G., Dere, K. P., Young, P. R., Landi, E., &: Mason, H. E. 2015, ArXiv e-prints 

Dere, K. P., Landi, E., Mason, H. E., Monsignori Eossi, B. C., & Young, P. R. 1997, A&AS, 
125, 149 

Eeldman, U., &: Doschek, G. A. 1978, A&A, 65, 215 
Eeldman, U., & Doschek, G. A. 1979, A&A, 79, 357 

Harper, G. M., Jordan, G., Judge, P. G., Robinson, R. D., Garpenter, K. G., & Brage, T. 1999, 
MNRAS, 303, L41 

Hibbert, A., Brage, T., &: Fleming, J. 2002, MNRAS, 333, 885 
Johnson, C. T., Kingston, A. E., & Dufton, P. L. 1986, MNRAS, 220, 155 
Keenan, P. P., et al. 2002, MNRAS, 337, 901 

Keenan, F. P., Crockett, P. J., Aggarwal, K. M., Jess, D. B., &: Mathioudakis, M. 2009, A&A, 
495, 359 

Landi, E., Del Zanna, G., Young, P. R., Dere, K. P., Mason, H. E., & Landini, M. 2006, ApJS, 
162, 261 

Liang, G. Y., Badnell, N. R., &: Zhao, G. 2012, A&A, 547, A87 

Olluri, K., Gudiksen, B. V., & Hansteen, V. H. 2013, ApJ, 767, 43 

Peter, H., et al. 2014, Science, 346, 1255726 

Straus, T., Fleck, B., & Andretta, V. 2015, ArXiv e-prints 

Tayal, S. S. 1999, Journal of Physics B Atomic Molecular Physics, 32, 5311 




7 
















0.45 



Log,o ( Density/cm 

Figure 4: The solid line shows the theoretical variation of the Oiv A1399.8/A1401.2 ratio as a 
function of density computed at log T = 5.15 using CHIANTI 8. The blue crosses show the 
measured ratios and densities with l-cr error bars. Numbers indicate the spectrum from which 
the measurements were made (Table [3|). 
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A. Sanity check 

To demonstrate that the low and high density limits of the Oiv A1399.8/A1401.2 ratio are 
reached in the coronal hole and flare kernel data-sets, respectively, we show in Figure [5] the 
A1399.8 and A1401.2 line profiles over-plotted in velocity space, with the A1401.2 profile scaled 
by the measured line ratio (Table [5|). 



Figure 5: A comparison of the A1399.8 (black) and A1401.2 (red) spectral line profiles from 
spectrum 8 (left) and spectrum 5 (right). The A1401.2 profile has been scaled by the measured 
ratio of the two lines, taken from Table [5l 


9 














